library(leaflet)
library(maptools)
library(rgeos)
library('rgdal')
alsace <- readShapePoly("data/inondation_shp/68_ATLAS_INONDATION_ZONE_INONDABLE",proj4string = CRS("+init=epsg:2154"))
alsace <- spTransform(alsace, "+init=epsg:4326")
leaflet() %>%addTiles()%>%
addPolygons(data=alsace)
set.seed(1)
df <- sp::SpatialPointsDataFrame(
cbind(
7.3895+rnorm(10,0,0.1), # lng
48.0661+rnorm(10,0,0.1) # lat
),
data.frame(type = rep("a",10))
)
# sp::coordinates(df) <- ~Long+Lat
sp::proj4string(df) <- sp::proj4string(alsace)
points=sp::over(df,alsace)
p=df[is.na(points$IDENT)==FALSE,]
p
## coordinates type
## 2 (7.407864, 48.10508) a
## 3 (7.305937, 48.00398) a
## 5 (7.422451, 48.17859) a
## 7 (7.438243, 48.06448) a
## 8 (7.463332, 48.16048) a
## 9 (7.447078, 48.14822) a
## 10 (7.358961, 48.12549) a
leaflet() %>%addTiles()%>%
addPolygons(data=alsace)%>%
addMarkers(data=p)%>%
addCircleMarkers(data=df,color="red")
dep_int <- readShapePoly("data/DEPARTEMENT/DEPARTEMENT",proj4string = CRS("+init=epsg:2154")) # ici nous avons du recharger le fond de carte departement avec le code 2154 qui signifie le format RGF93 / Lambert-93 de notre fond de carte de base, RGF93 :Réseau géodésique français légal pour la métropole depuis le décret du 26.12.2000
dep_int <- spTransform(dep_int, "+init=epsg:4326")
leaflet() %>%addTiles()%>%
addPolygons(data=dep_int)
sp::over(df,dep_int)
## ID_GEOFLA CODE_DEPT NOM_DEPT CODE_CHF NOM_CHF
## 1 DEPARTEM0000000000000005 68 HAUT-RHIN 066 COLMAR
## 2 DEPARTEM0000000000000005 68 HAUT-RHIN 066 COLMAR
## 3 DEPARTEM0000000000000005 68 HAUT-RHIN 066 COLMAR
## 4 DEPARTEM0000000000000005 68 HAUT-RHIN 066 COLMAR
## 5 DEPARTEM0000000000000005 68 HAUT-RHIN 066 COLMAR
## 6 DEPARTEM0000000000000005 68 HAUT-RHIN 066 COLMAR
## 7 DEPARTEM0000000000000005 68 HAUT-RHIN 066 COLMAR
## 8 DEPARTEM0000000000000005 68 HAUT-RHIN 066 COLMAR
## 9 DEPARTEM0000000000000005 68 HAUT-RHIN 066 COLMAR
## 10 DEPARTEM0000000000000005 68 HAUT-RHIN 066 COLMAR
## X_CHF_LIEU Y_CHF_LIEU X_CENTROID Y_CENTROID CODE_REG
## 1 1024125 6784581 1016174 6763894 44
## 2 1024125 6784581 1016174 6763894 44
## 3 1024125 6784581 1016174 6763894 44
## 4 1024125 6784581 1016174 6763894 44
## 5 1024125 6784581 1016174 6763894 44
## 6 1024125 6784581 1016174 6763894 44
## 7 1024125 6784581 1016174 6763894 44
## 8 1024125 6784581 1016174 6763894 44
## 9 1024125 6784581 1016174 6763894 44
## 10 1024125 6784581 1016174 6763894 44
## NOM_REG
## 1 ALSACE-CHAMPAGNE-ARDENNE-LORRAINE
## 2 ALSACE-CHAMPAGNE-ARDENNE-LORRAINE
## 3 ALSACE-CHAMPAGNE-ARDENNE-LORRAINE
## 4 ALSACE-CHAMPAGNE-ARDENNE-LORRAINE
## 5 ALSACE-CHAMPAGNE-ARDENNE-LORRAINE
## 6 ALSACE-CHAMPAGNE-ARDENNE-LORRAINE
## 7 ALSACE-CHAMPAGNE-ARDENNE-LORRAINE
## 8 ALSACE-CHAMPAGNE-ARDENNE-LORRAINE
## 9 ALSACE-CHAMPAGNE-ARDENNE-LORRAINE
## 10 ALSACE-CHAMPAGNE-ARDENNE-LORRAINE
Copyright © 2016 Kezhan SHI